An exact model of conformal quintessence 

Luis P. Chimento*, 1 Alejandro S. JakubiV and Diego Pavon* 2 

1 Departamento de Fisica, Universidad de Buenos Aires, 1428 Buenos Aires, Argentina 
2 Departamento de Fisica, Facultad de Ciencias, 
Universidad Autonoma de Barcelona, 
08193, Bellaterra (Barcelona), Spain 

Abstract 

A non-minimally coupled quintessence model is investigated and the conditions for a stationary 
solution to the coincidence problem are obtained. For a conformally coupled scalar field and 
dissipative matter, a general solution possessing late acceleration is found. It fits rather well the 
high redshift supernovae data and gives a good prediction of the age of the Universe. Likewise, the 
cold dark matter component dominates the cosmological perturbations at late times albeit they 
decrease with expansion. 
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I. INTRODUCTION 



After some years of research, the accelerated expansion of the Universe appears to 
have gained further ground 111 but the nature of dark energy -the substratum behind this 

nn 

acceleration- remains as elusive as ever While one may expect that one or other of the 

n n 

dark energy candidates (a very small cosmological constant [4|, quintessence |5[, Chaplygin 
[as (J, tachyon field interacting quintessence 0,0], non- minimally coupled quintessence 
lOL etc) will finally emerge as the successful model, at the time being none of them is 
in position to claim such status. 

In introducing dark energy as a novel component most of these candidates encounter the 
so-called "coincidence problem," namely, "why are the energy densities of both components 
(dark energy and dark matter) of the precisely same order today?" As shown by the au- 
thors, this problem has a dynamical solution provided that the dark matter component is 



aynamical solution proviaea tnat tne d 
or interacts with the dark energy g| 



In such a case, it can be 



assumed to be dissipative 
demonstrated that the equations governing the cosmic evolution imply the stationary con- 
dition (Eq. (J1(J|) below) and that the system is attracted to a stationary and stable solution 
characterized by the constancy of both density parameters, i.e., fl m and fl^ tend to constant 
values at late times. 

The aim of this paper is to provide an exact quintessence model, non-minimally coupled 
to the Ricci curvature. Non-minimal coupling naturally arises in generalizing the Klein- 
Gordon equation from Minkowski space to a curved space -for a recent review on this subject 
and further motivations, see Ref. 13|- We believe it is rather reasonable to explore whether 
a non-minimal coupling of the scalar field acting as dark energy to the Ricci curvature may 
be of help to understand the present stage of accelerated expansion and shed some light 
into the nature of the dark component. As a first step toward weighing the contribution 
of the coupling to the evolution of the Universe we shall consider the simplest case of a 
non-minimally coupled quintessence with a constant potential. As it turns out, there is a 
stationary solution that ensures that for late times the ratio between both energy densities 
remains a constant while the Universe asymptotically approaches a de Sitter expansion. In 
this late regime, the density perturbations of large wavelength decrease and are dark matter 
dominated. Further, our results are consistent with the high redshift supernovae data and 
provide a very reasonable estimate of the age of the Universe. 
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Section II presents our model which, aside from the non-minimally coupled scalar field, 
introduces a dissipative pressure in the dark matter fluid; this pressure turns out to be key 
for the solution of the coincidence problem. Section III derives the statefinder parameters as 
well as the age of the Universe. Section IV studies the observational constraints on our model 
imposed by the high redshift supernova data. Section V investigates the long wavelength 
density perturbations. Finally, section VI summarizes our findings. As usual, a subindex 
zero indicates that the corresponding quantity must be evaluated at the present time. We 
have chosen units so that c = 87rG = 1. 

II. CONFORMALLY COUPLED SCALAR FIELD 

Let us consider a Friedmann-Lemaitre-Robertson- Walker (FLRW) spacetime filled with 
two components, namely, dissipative matter and a non-minimally coupled quintessence field. 
The equation of state of the first component is of baryotropic type p m = (j m — l)p m , 
where the baryotropic index of matter is restricted to the range 1 < 7 m < 2. In addition 
to this equilibrium pressure the matter component is assumed to have a non-equilibrium 
(dissipative pressure) n connected to entropy production. It should be noted that barring 
superfluids (as Helium superfluid), this quantity is ever-present in every matter fluid found 
in Nature and is negative for expanding fluids |Q. In the case at hand, it may either 
come from interactions within the dark matter, or the decay of dark matter particles into 
dark particles Q], or from the non-linear growth of cosmic structures Q|, and it proves 
crucial to solve the coincidence problem. Likewise, the equation of state of the quintessence 
component can be written as p^ = (7^ — l)p^ with 7^ < 1. 

The Friedmann equation and the conservation equations for the matter fluid and 
quintessence read 



3H 2 + 3^ = p m + P4> (if = 1,0,-1) 



(1) 




(2) 



3 



+ 3iJ0 + £R0 + ^P = O, (3) 



ng const a 

0, Q 



where £ denotes the non-minimal coupling constant and the Ricci curvature scalar R is 
related to the quintessence scalar field by 



[1 - £(1 - 6£)0 2 ] R=-(l- 6£)0 2 + 4V - G^^P- + (4 - 3 7m )p m - 3vr . (4) 



Likewise, the energy density and pressure of the quintessence field are 

p = ±0 2 + V{(j>) + 3£#0(#0 + 20), (5) 



±</> 2 -V(0)-£ 4H(j)<i) + 2<j) 2 + 200+ (2H + 3H 



2^2 



(6) 



Obviously, these two reduce to their minimally coupled expressions for vanishing £. In terms 
of the density parameters Q m = p m /(3H 2 ), = p^/iSH 2 ) and Qk = ~K/(aH) 2 , the set 
of equations (HJ)-© become 



f2 m + fL, + — 1, 



(7) 



(lm + SH (^02 +lm + y)n m = 0, (8) 



^ + 3# [|^+7^)^ = 0, (9) 



respectively. 
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The simplest solution to the system of equations (JSJ)- that solves the coincidence 
problem is that fl m = Q m o and = fl^o at late times, with fl mQ and Q^o constants. This 



automatically implies the stationary condition 
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7T 2H 
lm + ^ = 1 * = 'MP' (10) 



It is readily seen that on the stationary solution, Eq.()10|). one has K = thereby we shall 
focus on spatially flat FLRW spacetimes hence forward. 

We are interested in obtaining a simplified, analytically integrable model that still retains 
the essentials of non-minimally coupled approaches and leads to an accelerated phase of 
expansion at late times. This may be accomplished by choosing for the coupling constant 
the conformal value £ = 1/6 and a constant value for the potential, V{<p) = Vq. As a 
consequence, Eqs. (JTJ), (|1J) and Q take a very simple form on the stationary solution 



3H 2 = (1 + r) 



I(0 + #0)2 + y o 



(11) 



R = A(l + r)V , (12) 



cu 2 



+ 3#0+— = 0, (13) 



where r = Qm/Q^ stands for the density ratio, and u 2 = (4/3)(l + t)Vq. In deriving (fT2^1 
we have made use of the fact that the Ricci scalar is R — 6 (if + 2H 2 ). It is interesting 
to note that, notwithstanding the potential being a constant, the gravitational interaction 
induces an effective mass given by m 2 s = u 2 /2 in the effective potential V c s = (ct> 2 2 /4) + V 
of the generalized Klein-Gordon equation. Thereby, in a loose sense, one might associate a 
particle -the "conformalon"- with the field <ft. 

In order to integrate this system of equations it is expedient to introduce the conformal 
time 7] and define a new field ip as 
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dt I X 

Tj = I — , f = <pa, 
a 



(14) 



then, Eqs. (|TT|) - (fT3j) become 



3( a ') 2 = (l + r) 



a" = —a 3 , i//' = . 



(15) 



where a prime indicates derivation with respect to 77. 

The general solution of this system of equations is given by 



vc sinhcut, 



(16) 



a \J c sinh cut 



where c and b are arbitrary integration constants, and the initial singularity has been fixed 
at t — 0. Combining Eq. (fT7)|) with the Friedmann equation on the stationary solution, 
3H 2 = (1 + r)p0, we get the following expression for the conformal quintessence energy 
density 

P(f> (t) = V coth 2 ujt. (18) 



In the late time accelerated regime, its equation of state 

pM = - ?* <i9) 

becomes that of a cosmological constant, and both density parameters fl^it) and fl m (t) = 
r£ltj,(t) asymptotically approach constant values. 
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Inserting Eq. (|16p into Eq. ((TH) it follows that 



0J\/C 



[1 - exp(-ut)) l/2 , 



(20) 



where F is the elliptic integral of the first kind. The conformal time is a growing monotonic 
function of t that behaves like y/i close to the singularity and has a finite upper bound 



17(00) 



2K(l 



UJ\/C 



(21) 



where K is the complete elliptic integral and 2X(l/v2) ~ 3.7. 

From Eq. (fTTjl it is apparent that depending on the choice of the constants several cases 
arise. If b ^ 0, then oc \j\ft for i — ► while if b = 0, the field has an extremum at the 
initial singularity. On the other hand, at late times the quintessence field evolves toward 
the minimum of its effective potential. 

The dissipative pressure 



7T 



-7m + 



3 cosh cut 



(22) 



follows from Eqs. (fTUj) and (fTT)|): and the evolution of the ratio n/p m , depicted in Fig. [TJ 
shows that the relative relevance of the dissipative pressure grows with expansion. From the 
last equation, when 7 m < 4/3, it is seen that there is a critical time t c given by 

cosh 2 ut c = (23) 



that separates the epoch with 7r > from the one with 7r < 0. From Eqs. (|23j) and (|16|) we 
see that it corresponds to the redshift 

where a = (sinhcuto) -2 is a characteristic parameter of the model with t denoting the 



1 



1, 



(24) 
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"age" of the Universe. As is well-known, the second law of thermodynamics requires 7r to 
be non-positive for expanding fluids 14]. This means that our model is valid for t > t c 
only ; when t < t c it describes a fictitious positive pressure contribution implying that any 
realistic evolution of the universe would yield a corrected value for the critical time t c > t c . 
Hence we would obtain a corrected age of the Universe t = t — t c + i c > t . This model 
may describe the growth of dissipative effects within dark matter as a consequence of the 
development of density inhomogeneities ^(|. Its characteristic growth time is given by the 
inverse of the conformalon mass, i.e., by Vq and r. In the late time accelerated regime, the 
transport equation for the dissipative pressure becomes n ~ —3(bH. So, we obtain for the 
dissipative coefficient (b ~ 7 m r[Vo/(3(l + r))] 1//2 which satisfies > as demanded by the 
second law. 



III. COSMOLOGICAL PARAMETERS 



The acceleration of the Universe is usually evaluated by the dimensionless deceleration 
parameter q = —a/(aH 2 ), where q < 0, q = 0, q > describes an accelerating, a linearly 
expanding (or contracting), and a decelerating universe, respectively. The present value go 
does not uniquely characterize the current accelerating phase thereby different dark energy 
models can lead to the same value. Useful additional information is encoded in the statefinder 
parameters r and s, defined as [l^ , |3| 



(25) 



aH*> 3(g -|; 



It is to be hoped that the pairs (r, q) and (s, q) will provide an accurate description of the 
present dynamics of the Universe and give us some insight into the nature of dark energy. 
This is only natural because r and s are directly connected to the third order term in Taylor's 
expansion of the scale factor around its present value [3] . In the case of an expansion given 
by Eq. (fT^j) these parameters are found to be 



a 2 a 2 4 — a 2 



s 



in terms of the adimensional ratio 



«-^ = m- ( 27 ) 

Hence there is a single functionally independent cosmological parameter, and the relation- 
ships 

2 r — 1 

r = 2 + q, s = , (28) 

y ' 3 2r - 5 ' y J 



between them hold (the dependence r(s) is depicted in Fig. |2J), so that the history of the 
deceleration parameter completely describes the evolution of this universe. Since q — > 1 
when t — > and q — > — 1 when t — > oo, this model describes a transition from a non- 
accelerated era to an accelerated era in the present Universe (see Fig. EJ). As the accelerated 
phase begins at a time t ac , where smh(ujt ac) — 1, we have that a ac — 



t ac = 1 cosh' 1 v 7 ^ = 4.31 v/ \ + a Gyr , (29) 
uj h 



where h indicates the current value of Hubble's constant in units of lOOkm/s/Mpc and 
the corresponding redshift is z ac = cr~ 1//4 — 1. For cold dark matter (j m = 1), the ratio 



cosh(u;t ac )/ cosh(a;t c ) = a/3/2 shows that at the commencement of this phase the dissipative 
pressure was already negative. On the other hand, in virtue of Eq. (I27jl . the age of the 
Universe can be expressed by 



t = - sinh' 1 = 4.31 v/1 , + g sinh" 1 Gyr, (30) 
uj \/a h Jo 



implying that a must be lower than unity if the Universe is to be accelerated at present. 
The time span since the critical time to the present is given by 

to — t c = — ( sinh -1 —= — sinh -1 —=) . (31) 

uj V VsJ 
The stationary condition (|l()jl shows that 7^ = 7 m at the critical time. Provided that by 

then the matter is cold, i.e., j m — 1 (in the next section we will verify the consistency of this 
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premise) we can assume a smooth extension of our model towards earlier times, as a cold dark 
matter dominated era for the purpose of obtaining a corrected age of the Universe. Indeed, 
using Eqs. (J2T)|) and (J27)l we find q(z c ) = 1/2 so that the deceleration parameter of 

this two-stage Universe is continuous. Hence, imposing the continuity of the matter energy 
density at the critical time, we have H c = uj and t c = 2/ (3a;). This yields the corrected age 

1 / 1 1 2\ 

t = - sinh" 1 — - sinh" 1 — + - . (32) 

As shown in the next section, this simple estimate produces a rather satisfactory result. 



IV. OBSERVATIONAL CONSTRAINTS 

It appears that supernovae of type la (SNela) may be used as standard candles. Prop- 
erly corrected, the difference in their apparent magnitudes is related to the cosmological 
parameters. Confrontation of cosmological models to recent observations of high redshift 



supernovae (z ^ 1) have f 
an accelerated expansion 



rown a good fit in reg ions of the parameter space compatible with 



20 



21 



22 



3, 



24, 



2511 . We note, however, that models like ACDM 



and QCDM usually require fine tuning to account for the observed ratio between dark energy 
and clustered matter, while our conformalon model, as well as QDDM/QIM models Q], Q, 
simultaneously provides a late accelerated expansion and solves the coincidence problem. 

Ignoring gravitational lensing effects, the predicted magnitude for an object at redshift z 
in a spatially flat homogeneous and isotropic universe is given by [26] 

m(z) = M + 51og£> L (z), (33) 



where M. is its Hubble radius free absolute magnitude and T>l is the luminosity distance in 
units of the Hubble radius, 

Di = (1 + z) f d/ 7/W < 34 > 

In virtue of Eq. (|27|) we obtain a representation in terms of the elliptic integral of the first 
kind 
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We have used the sample of 38 high redshift (0.18 < z < 0.83) supernovae of Ref. 2l| . 
supplemented with 16 low redshift (z < 0.1) supernovae from the Calan/Tololo Supernova 



Survey [27]. This is described as the "primary fit" or fit C in Ref. 2l|, where, for each 
supernova, its redshift Zi, the corrected magnitude m ; and its dispersion <7j were computed. 
We have determined the optimum fit of the conformalon model by minimizing a x 2 function 

X ~2 > ( 36 ) 

i=i * 

where A^ = 54 for this data set. The most likely values of these parameters are found to 
be (cr, Ai) = (0.2041,23.93), yielding xLJ n df = 1-104 (N DF = 52), and a goodness-of-fit 
P(x 2 > Xmin) = 0.282. These figures show that the fit of conformalon cosmology to this 
data set is even better than the fit of the ACDM or QDDM/QIM models in spite of the fact 
that we have at our disposal just one free parameter, namely, a. 

We estimate the probability density distribution of the parameters by evaluating the 
normalized likelihood 

«°< M ^ (37) 

Then we obtain the probability density distribution for a marginalizing p(a,Ai) over Ai. 
This probability density distribution p(cr) is shown in Fig. E]and it yields a = 0.224 ±0.054. 
We next use Eq. (jSOJ) to obtain H t = 0.833 ± 0.045. Likewise, taking h = 0.7 ± 0.07 (cf. 
j3|) we get from Eq. (p!Tj) a period since the critical time to — t c — 7.46 ± 1.03 Gyr, hence 
a corrected age of the Universe from Eq. (|32|) of to — 12.7 ± 1.4 Gyr. This one standard 
deviation range for the corrected age of the Universe falls within the 95% confidence age 
range 11.2 — 20 Gyr derived from the age of the oldest globular clusters and it is f ully 
consistent with the recent estimation of 13.4 ± 0.3 Gyr reported by the WMAP team [29j 
though the latter was reached on the basis of the standard ACDM model. 

By resorting to Eq. flU)) we obtain q = -0.809 ± 0.043, f = 1.191 ± 0.043 and s = 
-0.109±0.030. Assuming j m = 1 we get from Eq. J21) the critical redshift z c = 0.931±0.123. 
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This figure is consistent with the assumption that matter is already cold at critical time and 
with the range of redshifts of the supernovae used in the fit. On the other hand, the 
accelerated expansion era begins at z nr = 0.467 ± 0.093. This value matches the estimation 
fr „,n ACDM and the two epoch .node, of Re, Q, where the deception paratneter is 
constant within each stage. 

From Eq. ([27)1 we obtain the frequency parameter to = 0.129±0.013 Gyr -1 , corresponding 
to a conformalon effective mass m e ff = 1.91 ±0.20 x 10~ 33 eV; and combined with the current 

3 yields V = 8.16 ± 1.69 x 10~ 3 Gyr~ 2 . Finally, from 



32 



density ratio r — 0.56 ±0.07 
Eq. (|2^|l we get for the current ratio 



it a 

Pr 



o 3(1 + a) 



-0.758 + 0.051, (38) 



implying that nowadays the dissipative pressure plays a rather prominent role. 

V. COSMOLOGICAL PERTURBATIONS 

This section considers the evolution of long-wavelength scalar perturbations of this model. 
We 8hal . fo.low the method employed by Perrotta and Baccigalupi in Q S based on 
the formalism developed by Hwang [3J| to describe the evolution of perturbations in the 
synchronous gauge. In this gauge the perturbed metric takes the form 

ds 2 = a 2 [-drj 2 + (5ij + h^dxUx*] , (39) 



where the tensor hij represents the metric perturbations and its Fourier transform can be 
written as 



fry (x, 77) = d ke 



]3;„ „ik x 



k i k j fr(k,r7) + (k i k j --^)6C(k,r/) 



(40) 



Here fr denotes the trace of the tensor fr y - and ( represents its traceless component -for the 
sake of brevity we will omit the arguments (k, 77) henceforth. 
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The perturbed Einstein equations read 



k 2 ( - -Hh! 
2 

h" + 2Hti - 2k\ 
h" + G(" + 2H(ti + 6C') - 2k 2 C 



a 2 8p 

T ' 

a 2 {p + p)6 
2 

— 3a 2 Sp , 

— 3a 2 (p + p)£ 



(41) 

(42) 
(43) 
(44) 



where 7i = a'/ a, and the perturbed density Sp, pressure Sp, velocity divergence 6 and shear 
£ take the form (for a detailed definition of these terms see 



velocity 

Q,Q) 



W i w 

OPm H 5 7^,0 _ 3 5— 

er 2 er 



2 cr 



6a 2 



(45) 



1 n ^ 5F" H5F> fp-p 2A; 2 \ F'h! 
-F* R5<1)+ — + — — + ( + — ] 5F 



Spm H n 1" n r , <!>'•• 



er er 



3a' 



9a 2 



(46) 



(p + p)6 



{Pm + Pm ) "n 



a?F 



8F' + H5F) 



(47) 



(p + p)S 



3a 2 F 



SF + 3^[C + 



(48) 



where F — 1 — (c/> 2 /6). There remains the perturbed Klein-Gordon equation: 



j» + 2 W + ( A; 2 + ia 2 i?") <ty = ^ + y «2 • 



(49) 



where the perturbed Ricci scalar reads 



It has been shown in Ref. 



SR = — (h" - 3Hti + 2k 2 () . 

. 3er v ' 



(50) 



38j | that the density contrast 5 = Sp/p at large scales grows as 



f] 2 ~ a during the matter dominated era previous to the critical time t c . Here we investigate 
the behavior of the large scale density perturbations in the asymptotically de Sitter era. In 
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this regime we have a ~ (c/2) 1//2 exp(a>£/2) and 



2 

i] — 77(00) = A77 ~ , (51) 



hence TC ~ — 1/Arj and <p ~ — (bu/2)Ar] for A77 — > _ , where we have made use of Eq. (JT7jl 
to obtain 0. The solution of the system of equations (}4~Tj) - (f5TH) at the lowest order in Arj and 
k 2 is readily found to be 



h~DAr), (^——Arj, 5p ~ £p m ~ — 5R ~ — — — A?/ , ~ ~ — — Ar/ (52) 
6 4 6 



50 ~ AiAr/ (53) 

where the integration constants Z) and Ai are functions of the wavenumber k. Thus we find 
that matter perturbations are dominant at large times and it holds for the perturbation of 
the energy density ratio 5r ~ 5 m = 5p m / p m . The density contrast decreases in the late time 
regime as 5 oc 1/a, so that it has a peak during the period when viscous pressure grows. 

VI. CONCLUDING REMARKS 

We have presented a model of late acceleration that fits extraordinarily well the high 
redshift supernovae data, yields a good prediction for the age of the Universe and solves the 
coincidence problem. The only free parameter in the supernovae fit takes a natural value, 
a = 0.224 ±0.054. 

The dissipative pressure it in dark matter is a key ingredient of our model, and it can 
attain comparatively large values. Such pressure may arise from the interaction of cold dark 
matter with itself or the annihilation and/or decay of this component. Different models of 
cold dark matter that may show these features have been proposed recently -see, e.g., the 



pedagogical short review of Ref. [39(. No doubt the dissipative pressure inherent to these 
models may have a profound cosmological impact j^ . 
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The quintessence scalar field has a constant potential and it is non-minimally conformally 
coupled to the Ricci curvature. In our view, the model possesses two appealing features, 
namely: (i) an exact and simple solution (Eqs. (jl6|) -(|20) ) h and (ii) notwithstanding the 
potential is a constant (and consequently plays the role of a cosmological constant), the 
nonstandard kinetic energy term and the dissipation in the matter allows a stationary regime 
where the ratio of the energy densities remains constant. Likewise we have calculated the 
statefinder parameters and have shown that on this regime they are functionally dependent 
so that in this case the deceleration parameter is enough to describe these solutions. 

The present model may describe the growth of dissipative effects within dark matter with 
some kind of selfinteraction as a consequence of the development of density inhomogeneities 
after a critical time when these inhomogeneities become large enough. The large scale 
cosmological density perturbations are seen to decrease in the asymptotic de Sitter phase, 
with the matter perturbations dominating over the quintessence perturbations. 

We believe that because of its simplicity, our conformal quintessence model may serve as a 
starting point for more complete models that include a nontrivial selfinteraction quintessence 
potential, consider the approach of the energy density ratio towards this stationary regime, 
and extend towards earlier times, when dissipative effects within dark matter become neg- 
ligible. 
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FIG. 1: Selected curves of the ratio dissipative pressure vs energy density of matter Tr/p m vs the 
redshift z between the present z = and z = 1. From top to bottom, the curves correspond to 
parameter a = 0.3, 0.224, and 0.18 defined by a = (sinhwto) -2 - 
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FIG. 2: Statefinder parameter r in terms of s, the other member of the pair. The evolution goes 
from the point (s c ,r c ) = (-co, 5/2), corresponding to the critical time to the point (sco^oo) = 
(0, 1), corresponding to the asymptotically exponential expansion at large times. The present state 
of the universe corresponds to (s ,f ) = (-0.109 ± 0.030, 1.191 ± 0.043). 
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FIG. 3: Selected curves of the deceleration parameter q vs the redshift z between the present z = 
and z = 1. From top to bottom, the curves correspond to a = 0.3, 0.224, and 0.18. 
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FIG. 4: The estimated probability density distribution (normalized likelihood) for the parameter 
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